DESY 01-161 
HU-EP-01/42 
October 2001 



1 
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We study the continuum limit of the abelian string tension and the density of abelian monopoles calculated 
after carefully fixing the maximal abelian gauge by employing the simulated annealing algorithm. We present the 
evidence that the abelian string tension converges to the nonabelian one in the continuum limit. For the monopole 
density we confirm earlier findings that the density of the properly defined infrared monopoles has correct scaling 
while the total density seems divergent in the continuum limit due to ultraviolate contributions. We also compare 
with results obtained with the usual iterative gauge fixing algorithm. 



1. INTRODUCTION 

Numerous lattice studies have provided sup- 
port for the dual superconductor scenario of con- 
finement ||l[. Most of the results have been ob- 
tained in the maximal abelian gauge (MAG) at 
particular values of the lattice spacing. To our 
knowledge no attempt was made for extrapola- 
tion into the continuum limit. In this paper we 
want to fill in this gap. 

2. GAUGE FIXING AND SIMULATION 
DETAILS 

We work with SU{2) gluodynamics in the MAG 
which is defined by maximizing the functional 

^(f^') = ^ E [t^.^^3C/.^>3] (1) 

with respect to the gauge transformations : 
Ui.f_, = 9xUx,tj,gl+^- The abehan projection 
Ux.^j, Ux.p, is defined by the following relations 

[/::^ = cos(0,,^)e»^i.., el^^-el^^ex,^. (2) 

Ux,f_i =diag {e'^- '', e"'^="-"}. (3) 

It is well known that the functional (Q) has many 
local maxima corresponding to Gribov copies . 
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0.12 fm 
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L ■ a 


2.9 fm 


2.1 fm 


1.7 fm 


< F >SA 


.7336(1) 


.7511(1) 


.7662(1) 


< F >Ro 


.7310(1) 


.7491(1) 


.7648(1) 



The standard procedure consists in selecting a lo- 
cal maximum one finds more or less randomly by 
an iterative algorithm. Our prescription of gauge 
fixing is to find the global maximum of (^. This 
prescription resembles the minimal Landau gauge 
§. The 'global maxima' definition of gauge fixing 
is known to be a limit of the global gauge fixing 
introduced in . As a numerical tool we use the 
simulated annealing (SA) algorithm This 
algorithm helps to get closer to the global max- 
ima but at some low auxiliary temperature T it, 
nevertheless, gets trapped into metastable states 
of the corresponding spin glass. To overcome this 
problem we applied the procedure to every MC 
configuration 10 times by starting from different 
gauge copies and used the copy with maximal 
value of F{U). For comparison we carried out 
gauge fixing with the standard iterative algorithm 
with alternating relaxation and over-relaxation 
sweeps (RO algorithm). Only one gauge copy was 
considered as it is commonly done. 
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The simulation details are presented in Table 1. 
There we also show the average value oi F{U) for 
the two algorithms. We investigated 20 statisti- 
cally independent configurations per /3-value. 

3. MONOPOLE DENSITY 

We use the DeGrand-Toussaint definition of 
magnetic currents |^ 



47r 



(4) 



with Ox^ufi = dij_9x^u - duOx.^i + 27rma;^p,y and 



-TT < 



.11/ < TI", 



0,±1,±2. 



The physical monopolc density is then defined as 
1 



(5) 



Following 1^ we will call monopoles relevant for 
confinement 'infrared' (IR) monopoles and the 
rest of monopoles, normally forming small loops, 
'ultraviolate' (UV) ones. The corresponding den- 
sities are denoted as pl^^^ ■ In jl^ it has 
been found that monopoles from the largest clus- 
ter (LC) are sufficient to reproduce the string 
tension. This observation suggests that the IR 
monopoles can be identified as the monopoles be- 
longing to the LC 1^. However, as it has been 
indicated in the LC may split into parts if 
the volume is not large enough. We confirm this 
observation. Moreover, we have found that, when 
the LC splits, then there usually appear clusters 
with nontrivial winding 



u^cluster 



(6) 



Therefore, we suggest the following definition of 
IR monopoles: monopoles from LC plus those 
from clusters with nontrivial w^. In Fig. |l| we 
present our results for pi^'^^, measured accord- 
ing to the new definition. ^ We also show our 
results for the density of monopoles belonging to 
the LC, /Onion, together with data from We 



^ Throughout this work we use_th£_following values of the 
nonabeUan string tension Vari? 0.266(2) at /3 = 2.4; 

0.1905(8) at 13 = 2.5; 0.136(4) at /3 = 2.6. We set the scale 
putting ^ = 440 Mev. 
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Figure 1. Density of monopoles vs. lattice spac- 
ing. 

found that when the lattice size is large enough 
the clusters with nontrivial windings are rare and 
the two definitions give coinciding results. This 
is also supported by the agreement between our 
value obtained at (3 = 2.5,L ■ a = 2.1fm 
and the Pmon value from obtained at = 
2.5115, L ■ a = 2.8 fm. We conclude that the new 
definition of IR monopoles can be successfully 
used on lattices with size L • a > 1.7 fm. It is to 
be checked, how this definition works on smaller 
lattices. We extrapolated to the continuum 
limit using a quadratic fit (solid curve in Fig. 1). 
We find that p^^„ = 0.65(2)cr3/2. Data in Fig. 1 
show that while is finite, p^on seems to di- 

verge providing a divergent total density. Our 
data indicate that p'. 
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. ABELIAN STATIC POTENTIAL AND 
STRING TENSION 

We introduce the abelian Wilson loop 



WabiC) = -Tr 



n 

\iec 



ui 



and the abelian static potential 

{WabiR. T)) 



{Wab{R,T+l)) 



(7) 



(8) 



On- axes as well as (0,1,1) off-axes directions were 
used. We apply smearing for spatial abelian 
links to improve the signal-to-noise ratio. Vj5J(i?) 
has been fitted to a constant within the range 
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Figure 2. Ratio of the abelian and nonabelian 
string tensions vs. lattice spacing for SA and RO 
algorithms. The solid line represents a quadratic 
fit. For comparison data from are shown. 



5 < Ta < 8 to obtain Vab{R)- We employ a di- 
mensionless form of the usual parametrization for 
the static potential 



Vab{R) 



RV^ 



a 



(9) 



Data points were fitted to this form in the range 
0.5 < i? • i/a < 2.5. Fig. || presents our main re- 
sult. The continuum limit of Uah agrees with that 
of a within error bars. One can see a good agree- 
ment with a result from obtained by the same 
method. It is worth noting that our result is still 
biased by Gribov copies effects, although these 
effects were substantially reduced in comparison 
with the common procedure. Using results of 
we can estimate the bias as t^Uahj^y ~ —0.03. 
Taking into account this bias would not spoil our 
conclusion. 

We also made an attempt to compute the 
charge-two abelian static potential. In |^ 12| it 
has been conjectured that the charge-two abelian 
string tension cr^^^ is equal to the adjoint string 
tension. If this is true, then in the continuum 
limit o'^^'^ ja 8/3. Unfortunately, our statistics 
appeared to be insufficient to provide a reliable 
estimation of c^jj"^ • The work on this and also on 
the monopole static potential is in progress. 



5. CONCLUSIONS 

- With Gribov copies effects taken under bet- 
ter control we found that in the continuum 
limit CTab/c — 1.03(4) supporting abelian 
dominance for the string tension within the 
maximal abelian gauge. 

- The density of IR monopoles turns out 
finite in the continuum limit: p^f^„ = 
0.65(2)cr3/2. Thus, the average distance 
between monopoles can be estimated as 
dmon = 1.15(1)/V^ ^ 0.52(1) fm. 
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